IFUP-TH 1/96 
January 1996 



lquant-ph/9601024 



The interaction time of a packet with a 
potential barrier 

Andrea Begliuomini" and Luciano Bracci a ' & 

(a) Dipartimento di Fisica dell'Universita, 
piazza Torricelli 2, 1-56126 Pisa, Italy. 
(b) Istituto Nazionale di Fisica Nucleare, sezione di Pisa, Italy. 



Abstract 

We study the evolution of a wave packet impinging onto a one dimensional 
potential barrier. The transmission and reflection times discussed in the lit- 
erature for stationary states do not correspond to the times required for the 
emergence of a transmitted or a reflected packet. We propose new definitions 
for the interaction (dwell) time and the transmission and reflection times which 
are suitable for packets and fit better the actual time evolution of the packet. 



PACS numbers: 03.65.Bz,73.40.Gk 



I. INTRODUCTION 



The tunneling of a particle beyond a potential barrier is one of the simplest effects 
predicted by Quantum Mechanics where the conflict between classical and quantum pic- 
tures is most striking. Although it has been studied since the early days of Quantum 
Mechanics |l| , the debate is still open as to what should be meant for tunneling time (see 



It could be objected that in a proper quantum formulation of the problem there is 
no room for such a concept, all that we can ask being the probability of detecting the 
particle beyond the barrier or of having it reflected by the potential. We feel, however, 
that this is an extreme view. It is legitimate, for example, to think of an ensemble of 
systems prepared in a given initial state, for which the times are measured when a detector 
located in front or beyond the barrier reveals the arrival of the particle (in order not to 
disturb the state, we think of separate experiments for measuring the arrival and the exit 
times). Knowledge of the initial state (i.e. the wave function at time t = 0) should enable 
us to determine which is the difference between the average exit time and the average 
arrival time. This difference could be regarded as the average tunneling time. While we 
do not advocate the above definition as the definition of the tunneling time, we insist that 
the wave function contains implicitly the information as to the time that in the average 
the particle spends in the potential barrier. The problem with the determination of this 
average time is that, the process being intrinsically non-classical, it is not possible to look 
for a quantum counterpart of a classical observable, whose average value on the state of 
the particle should be interpreted as the tunneling time. On the other hand, a clearer view 
about the tunneling time is urged also by current experiments on semi-conductor devices 
H and evanescent waves 0], where observation of superluminal velocity has sometimes 
been claimed. 

Most of the approaches to the problem of tunneling time deal with stationary states. 
The particle is in an eigenstate of the Hamiltonian, and the tunneling time tt, as well as 
the reflection time tr and the dwell time td, are functions of its energy E, or momentum 



k = ^2mE/h. Also those who envisage wave packets (see [§], 0, ||) use packets which 

are so narrow in energy as to allow the monochromatic approximation to hold. We have 
preferred to investigate the time evolution of a wave packet with a spread in energy which 
forbids this approximation. The packet impinges onto the barrier and is partially reflected 
and partially transmitted. It is constructed as a Gaussian superposition of eigenstates 
centered around a value E lower than the height Vq of the barrier. We observe the 
evolution of the packet at times tt(E), Tr(E) and Td(E) after it has reached the barrier. 
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Although the time t in when the packet "reaches the barrier" is not sharply defined, due 
to the interference of the higher energy components, which reach the barrier earlier and 
are partially reflected, with the lower energy incoming components, we find that the 
above mentioned times tt(E), t&(E) and td(E) definitely do not correspond to the times 
required for the packet to emerge from the barrier. The uncertainty in t m is not such as 
to alter this conclusion. It can be argued that the transmitted packet is centered around 
an energy Et higher than E, whereas the reflected packet is centered around an energy 
£"r lower than E, but also the times tt(-Et) o r t r{Er) are by no means representative of 
the times required to see the reflected packet or the transmitted packet. The conclusion 
seems to be that the tunneling times defined for stationary states with energy E are not 
meaningful for the evolution of a wave packet having E as average energy. 

This has prompted us to look for another determination of the times tt, 7d and tr 
which take into account the actual behaviour of the wave packet. As for td, we define 
it as the time integral of the probability P 2 (t) of finding the particle within the barrier 
region, 

TD = / P 2 (t)dt (1) 

this definition being unsensitive to the time t 2 when the particle begins to interact with 
the barrier, provided it is earlier than the time when the packet impinges onto the barrier. 
Thus, Td can be regarded as the total interaction time. As for tt, we define it as a weighted 
sum of every time interval At, the weight being the fraction of the transmitted packet 
which at time t has not yet been transmitted. The definition for tr is similar. The trouble 
with these definitions is that there is still a problem with the time when the integral over 
t begins. In principle, the packet is interacting with the barrier since time t = 0, but the 
contribution of the earlier times, when the particle has not yet arrived onto the barrier, 
to r T should be negligible. 

While there is no objective way of determining the time when the interaction between 
the particle and the potential starts, we observe that a shifting of t 2 in eq. (|I|) affects 
td, whose value decreases with increasing t 2 . On the other hand any determination of td 
entails a given accuracy e. For the lower integration limit in the definitions of tt and tr 
we choose the value ti, such that // 2 L P 2 (t)dt equals e. The meaning of this choice is that we 
have a finite time resolution e, and we neglect those time intervals which contribute to the 
interaction time less than e. We note however that the time dependence of the probability 
P 2 (t) and of the analogous probabilities of finding the particle beyond the barrier (P 3 (t)) 
and in front of it (Pi(t)) when the packet impinges onto the barrier is sufficiently steep 
so that the choice of £l does not really affect tt and tr in a substantial way. The values 
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we find (with e ~ 0.01) for tt, tr and Td are definitely different from the values of tt(E), 
tr(E) and td(-E') and are in agreement with the actual time evolution of the packet. 

We have also examined the behaviour of the probabilities Pi{t) for large t. These 
probabilities have an exponential tail e _<//r , with the same time constant r (depletion 
time) for each probability. The value of r is determined by the position of the poles of the 
transmission coefficient D(k) in the complex k plane, and is independent of the details of 
the wave packet. 

In conclusion, we find that the tunneling times proposed in the literature for stationary 
states are not meaningful for the actual time evolution of a wave packet. The time lapse 
from the first contact of the packet with the barrier to its emergence beyond the barrier, 
although not so sharply definable, is definitely different from the tunneling times proposed 
for stationary states. 

In section 2 we define the problem and build the packet whose time evolution is 
discussed in section 3. In section 4 we present the definitions of td, tt and tr and in 
section 5 we discuss the depletion time. Section 6 is devoted to the conclusions. 



II. THE WAVE PACKET 



We consider a one dimensional problem, the Hamiltonian being 



*=£ + V(,) (2) 



. V (\x\ < d) 

v(x) = < o u > d) (3) 

We use units such that H = 1. The eigenf unctions ipk of the Hamiltonian (fj) are well 
known (see Appendix). We consider the evolution of a wave packet ip built as a Gaussian 
superposition of the functions ip^i impinging onto the barrier: 

V>(ar, t) = J a{k) ^ k {x) e - ikH / 2m dk (4) 

where the coefficients a(k) are 

k av is the average momentum, and xq is the coordinate of the peak of the packet at t — 0. 
We choose the parameters involved in the problem as follows: 



3 



m = \ K v = 9.9 5 = V2 x = -15 d = 2 K = 50 (A*, = 10) (6) 

The packet flD is Gaussian also in x (fig. 1). The reason why the peak of the packet is 
located so far from the left edge of the barrier at time t = is to have an identification 
of ti n as sharp as possible. With a packet starting nearer to the barrier, the Gaussian 
form of the packet would be immediately lost, due to the interference of the incoming and 
reflected components of the packet. 

III. EVOLUTION OF THE PACKET 

We have studied the evolution of the packet in order to verify to what extent the 
definitions of the tunneling time proposed for stationary problems are meaningful for 
a wave packet. More precisely, we have tested whether the phase times r ph |7j or the 
times proposed by Buttiker r B || (which turn out to be deeply connected with the times 
proposed within other approaches, see for ||- 0] a review) do represent the lapse of time 
which the packet spends in the barrier. To this purpose, it is necessary to mark the time t; n 
when the packet begins to interact with the barrier, and the time t nn when a transmitted 
(reflected) packet appears. The comparison of t nn — t in with r ph and r B calculated for 
significant values of the energy will show that r ph and r B are not significant for a wave 
packet. 

We observe the shape of the packet as it moves towards the potential barrier. As 
long as it is far enough, its shape is quite similar to the initial shape; when it approaches 
the left edge of the barrier it begins to become blurred (fig. 2) due to the interference 
between the incoming and the reflected components. There is a time interval in which 
the interference phenomenon is dominant, but a reflected packet is still absent. The peak 
of a reflected packet appears a time At after the blurring of the incoming packet. A bit 
later, we see the emergence of a transmitted packet beyond the barrier. 

For the time t; n when the packet begins to interact with the barrier we assume the 
time when the blurred shape can be macroscopically observed. As for the time t nn , we 
have two possibilities: the time when the reflected packet appears (t n n,R,) or the time when 
the peak of the transmitted packet appears (tfin/r)- ^ is clear from the above that neither 
tin nor t nn ,R (£fin,T) are sharply defined. In the units we have chosen, each of them can be 
determined only within an error At ~ 0.1. 

By inspection of the graphs representing the reflected and the transmitted packet 
respectively (see figs. 2a and fig2b), we get 

tin ^0.9 £fin,R-l-9 t fin ,T = 2.7 (7) 
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For the times tr and tt it follows 

r R ~ 1, r T ~ 1.8 (8) 

The error on the above values can be assessed to be of order 0.2. 

We compare the above times with the times r ph and r B . These are functions of the 
momentum k, so we must decide which momentum to consider. We consider the average 
momentum k av and the average momenta &r and fcx of the reflected and transmitted 
components respectively. These latter are calculated to be 

k R = 9.696, k T = 10.327 (9) 

Incidentally, this shows, as previously noted that a potential barrier acts as an accel- 
erator: the transmitted wave packet has an average momentum larger than the incoming 
one. The opposite holds for the reflected packet. 

In table 1 we present the dwell times t ph ' D and t B,D calculated for k = k av , &r and k^, 
together with the reflection time t B,R calculated for k av and &r and the transmission time 
t B,T calculated for k av and k^. We see that the dwell times t ph > D and t B,D are definitely 
shorter than the tunneling time tt reported in (|8|). The same holds for the reflection 
times t B,R calculated for /c av and &r. As for the transmission time t B,T , we see that the 
value corresponding to k av is longer, whereas the value corresponding to k^ is shorter. We 
conclude that the tunneling times found for stationary problems are not useful to describe 
the evolution of a packet. 



IV. DEFINITION OF r D , r T AND r R 

So far definitions of the time that the particle interacts with the potential (the so called 
dwell time tq) and of the transmission (tt) and reflection (tr) times have been given 
mainly for stationary problems. Even the authors who have dealt with wave packets 
considered packets which were so narrow in energy that the relevant times could be 
considered to be functions r(k) of the momentum, as in the stationary case. In this 
section we propose definitions of td, tt and tr which are suitable for a wave packet. 

We first consider the probability Pi(i) that at time t the particle is within the barrier 
region, 

P 2 (t)= I* dx \tP(x,t)\ 2 (10) 
J-d 

and the analogous probabilities that at time t the particle is in front of the barrier (Pi(t)) 
or beyond the barrier (P 3 (t)). Obviously, we have P\{t) + Pzif) + Pzif) = 1 
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The values of Pi(t), P%{t) and Ps(t) are reported in fig. 4. P x and P 3 tend to asymptotic 
values which we call respectively R and T. They represent the probabilities that the 
particle is reflected or transmitted respectively. Obviously, R + T = 1. 

Inspection of P2{t) shows that for £ < 0.75 the packet does not interact with the barrier, 
and the interaction reaches its maximum at £ ~ 1.5. After this time the probability of 
finding the particle in the barrier region decreases with a tail which has an exponential 
shape. The interaction time should be the total time that the particle spends in the 
potential region. With this view, we propose the following (already deined as eq. ([!])) as 
a definition of Tp (see also ref. [TH], [|TT|] where a similar definition is put forth): 

r+oo 

td = / P2{t)dt 

The meaning of eq. ([!]) is clear: every time interval dt is weighted with the proba- 
bility p2(t) of finding the particle within the potential barrier. The dwell time is to be 
interpreted as the time the particle interacts with the potential regardless its fate. The 
definition is independent of the choice of time t2, provided it is chosen earlier than the 
time the packet is significantly present in the barrier region. We can safely take £2 = 0. 

We define the transmission time tt as the average time that it takes for the transmitted 
particles to emerge beyond the potential barrier. This time is given by the integral 

(1 - Z®) * (id 

In the above definition 1 — P$(t)/T is the fraction of the transmitted packet which has 
not yet been transmitted at time £. We interpret this fact viewing this fraction as "being 
transmitted" at time £, so that any time interval dt contributes to the transmission time 
with a weight 1 — P 3 (£)/T. 

The trouble with this definition is the lower integration limit £3, in that eq. ( |i~T|) 
gives weight 1 also to the time intervals when the packet has not yet arrived in the 
potential region. But these time intervals are not to be considered as contributing to the 
transmission time: what eq. (|ll]) actually gives is the transmission time starting from 
time £3. 

The problem cannot be circumvented by any choice of time £3. In principle, £3 should 
be the time when the particle begins to interact with the potential, but this time cannot 
be determined in any objective way. However, we can determine that time t e such that 
the packet has spent in the potential region an amount e of time starting from the time 
£ = 0. This is the time t e such that the integral of P2 from £ = to £ = t e equals e. Now, 
the transmission time Vr(e) reckoned from time t e (such that the time spent in the barrier 
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by the packet is equal to e) is denned unambiguously. On the other hand, the calculation 
of the dwell time td (as well as any possible time measurement about the particle) is 
affected by an error. If we choose e to be the same as this error, the uncertainty on r T 
due to the choice of the lower integration limit can be considered of the same order as the 
time resolution we are able to attain. In conclusion, we put 



/OO / ft 

6 j P 2 (x)dx-e 



1 - 



T 



dt 



(12) 



where @(x) is the Heavyside step function. The lower integration limit in eq. (|12"|) can be 
taken the same as in eq. (|TT|). 

Along these lines we define also the reflection time Tj^(e); 



r+infty / rt 

/ 6f / P 2 (x)dx-e 



1 - 



R 



dt 



(13) 



The interpretation of eq. ( |1~3"D is straightforward. 1 — P\{t)/R is the fraction of the 
reflected packet that at time t has not yet been reflected. This fraction is taken as the 
weight for any time interval dt. We note however that for times near the beginning of the 
interaction of the particle with the potential this weight can be negative. But the decrease 
of Pi{t) is very sharp (see fig. 4) and the contribution to the integral of the region where 
the weight is negative is small indeed. 

In the case we have investigated we have found 



r D = 0.93 



In order to find r T and r R we have fixed e = 10 r D ~ 0.01. This yields 



r T = 3.39 tr = 0.55 



(14) 



(15) 



Note that we have T ~ 0.14, R ~ 0.86. With this values and the results reported in eqs. 
14 and 15 the conditional probabilty relation ( ||, [[11 1 



td = Tr T + Rt] 



(16) 



is fairly satisfied. Condition (|T6| ) would be identically satisfied if the lower limits where 
the integrands in eqs. (H), ( |i~2"D and ([13]) start to differ from zero were the same. The 
fact that eq. (|16|) holds true with a fair accuracy, to within that value .01 which can 
be assessed as the accuracy of all our calculations, can be regarded as a support to the 
correctness of the definitions of td, tt and tr. 



7 



By inspecting table 1 we see that td is definitely larger than r p h,D and t Bi d evaluated 
for k = k av . The dwell time td looks a very reliable estimate of the interaction time. The 
discrepancy with the value found for stationary states confirms that the extrapolation to 
wave packet is untenable. As for tt and tr, a comparison with the previously reported 
values of the reflection and transmission times derived by inspection of the wave packet 
evolution shows that r T and r R in eq. (^) are respectively longer and shorter. Thus, the 
discrepancy between the times in eq. ( |I~5D and the times derived within the stationary 
approach, for example the Buttiker times calculated for k = k av , is even larger. 

V. THE DEPLETION RATE 

The tail of P2{t) can be described as an exponential curve with a time constant Td cp : 

P 2 (t) = Ae~ t/Td ^ for large t (17) 

By considering the values of P2 for t > 30, we find that the exponential fit is excellent, 
with 

r dep = 16.192 (18) 

and a correlation coefficient R = —0.9999883. The same time constant Td ep rules the 
asymptotic behaviour of Ps(t) and, due to probability conservation, of Pi. 

r dep = 16.192 R = -0.9999435 

The value of T dep is connected with the behaviour of the complex transmission coefficient 
D(k) in the complex plane. By explicitly writing ( see Appendix) we see that the 

only singularities are in the product u(k)u*(p) in the denominator, with 

u(k) = (k 2 - k 2 ) sinh(2/«f) - 2ikn cosh(2/cd) (19) 

It is easy to see that if u(k) = 0, then u*(k*) = 0. Hence, a zero for u in k = x + iy, 
together with the zero in x — iy for u*, will contribute to Pi(t) with a term exp(2xyt/m). 
For large values of t, the term with xy negative and minimum in absolute value will 
dominate. The behaviour of P 2 will be as in eq. (|TTD, with: 

T dep = m/2\xy\ (20) 

The search for the pole with smallest \xy\ value is discussed in the appendix. We find 
x = 10.03, y = —3.0565 • 10~ 3 , which yields r dep = 16.3087, in fair agreement with the 



value reported in eq. (|T8D . 
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VI. CONCLUSIONS 

We have studied the evolution of a wave packet which approaches a potential barrier. 
The time requested for the appearance of a transmitted (tt) or a reflected (tr) packet, 
reckoned from the moment the incoming packet begins to interact with the barrier, are 
definitely different from the values of tt and tr reported in the literature for monochro- 
matic packets. We propose new definitions of td, tt and tr which are suitable for a packet 
and give results that fit better the actual time evolution of the packet. 

The probability of finding the particle within the barrier has an exponential tail whose 
time constant Td ep is determined by the behaviour of the stationary solutions in the 
complex momentum plane. 



APPENDIX: 



The stationary solutions of the Schroedinger equation with the Hamiltonian (g) are: 

e ikx + A(k)e~ ikx if x < -d 
V> fc (x) = { B(k)e KX + C(k)e~ KX if \x\ < d 
D(k)e ikx if x > d 



where 



A ^ = (/t 2 + A; 2 )sinh(2^) c _ 2tfcd 
u(k) 

= _ ik(K + ik) e _ iMe _ Kd 

u{k) 

= Jk{K - ik) e _ ikdeKd 

u(k) 

D(k) = -^ e ~^ 
u{k) 



k = ^k 2 Q - k 2 
u(k) = (k 2 — k 2 ) sinh(2K(i) — 2ikn cosh(2K(i) 
The probability is given by: 



p 2 (t) = JP[dk[ dpe-( p - k ^ 2s2 e-( k - k ™) 2 s 2 e -^ 2 -p 2 Wi™ 



2tt 3 



kpe 



i(p—k)d 



[(q 2 — p 2 ) sinh(2g<i) — 2ipq cosh(2qd)]* {(k 2 — k 2 ) sinh(2/«c?) — 2ikn cosh(2/?c 



( 1 _ e -2(K+q)d\ ( 1 _ e -2(K-q)d\ 

(q + ip)*(K + ik)± = + (q- ip)*{K + ik)± = J -+ 



— (q + ip)*(K — ik) 



k + q 

_ e 2(K-q)d 

k — q 



— (q — ip)*(n — ik) 



k — q 

K + q 



where q — y k% — p 2 

The only singularities are in the factors u(k)u*(p) in the denominator. In order to find 
the zeroes of u(k) we note that if u(k) = 0, then 



sin 



2ik VT^ r z 2 ~ 



4z 2 



(Al) 



with z = k/ko. For an opaque barrier (i.e. 2k d ^> 1) the solution for z has to be near 
z — 1. The values of the zeroes have been found by solving for z the relation 

= u(z) ~ u (z ) + (z - Zq) u (z ) 

and iterating, choosing different values of z near 1. We have found as many different 
solutions as predicted by the Cernlib routine nzeros for a neighbourhood of z — 1. The 
value of Td e p is given by that solution for which xy is negative and minimum in absolute 
value. 
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FIGURE CAPTIONS 

Figure 1 

The Gaussian shape of the incoming packet for t — 0. 



Figure 2a and 2b 

In (a) the wave packet at t — 0.9 and in (b) a zoom of the blurred region. 



Figure 3a and 3b 

In (a) the reflected packet that appears at t — 1.9. In (b) the first transmitted peak that 
appears at t — 2.7. 



Figure 4 

Plots of Pi(t), i"2(0 and Ps(t) as functions of t. 



TABLE CAPTIONS 

Table 1 

The phase time and the Buttiker times evaluated in k — k, kn, kr- 
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time 


k — k &v 


k = k R 


k = kx 


r ph,D 


0.143 


0.0843 


1.011 


r B,D 


0.140 


0.079 


1.008 


T B,T 


2.357 




1.248 


T R 


0.140 


0.079 





Table 1 
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